Gauss jacobi method